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Abstract 

We present an approximate, analytical calculation of the reionized 
spectra Cf^^ of cosmic microwave background radiation (CMB) 
anisotropics and polarizations generated by relic gravitational waves 
(RGWs). Three simple models of reionization are explored, whose 
visibility functions are fitted by gaussian type of functions as ap- 
proximations. We have derived the analytical polarization /?/ and 
temperature anisotropics ai, both consisting of two terms propor- 
tional to RGWs at the decoupling and at the reionization as well. 
The explicit dependence of Pi and ai upon the reionization time 
ijr: the duration Ar/^, and the optical depth Kr are demonstrated. 
Moreover, Pi and ai contain Kr in different coefficients, and the po- 
larization spectra C^^ are C^^ are more sensitive probes of reion- 
ization than Ci^^. These results facilitate examination of the reion- 
ization effects, in particular, the degeneracies of Kr with the nor- 
malization amplitude and with the initial spectral index of RGWs. 

*yzh@ustc. edu.cn 
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It is also found that reionization also causes a /t^-dependent shift 
Al ~ 20 of the zero multipole Iq of Cf^, an effect that should be 
included in order to detect the traces of RGWs. Compared with 
numerical results, the analytical C/^^ as approximation have the 
limitation. For the primary peaks in the range I ~ (30,600), the 
error is < 3% in three models. In the range / < 20 for the reion- 
ization bumps, the error is < 15% for C^^ and C^^ in the two 
extended reionization models, and C/^-^ and Cf^^ have much larger 
departures for / < 10. The bumps in the sudden reionization model 
are too low. 

PACS numbers: 98.70.Vc, 98.80.-k, 04.30.-w, 04.30. Nk, 
Key words: 
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1. Introduction 



Reionization is a very important cosmological process, which might be, to a large extent, 
determined by the first luminous objects formed in the early universe, either star-forming galaxies 
or active galactic nuclei. Our knowledge of the cosmic structure formation of the universe would 
be incomplete without a reliable account of reionization history, the details of which is still not 
understood yet. During the evolution history of CMB, the reionization taking place around the 
redshift z = (20 ~ 6) is a major process in shaping the profiles of CMB spectra on large scales, 
only secondary to the decoupling around z ~ 1100. Reionization leaves observable prints on CMB 
[1, 2, 3, 4, 5, 6, 7, 8] through the interaction between the CMB photons and the reionized free 
electrons. In particular, the spectra of CMB anisotropies and polarizations on large angular scales 
contain the distinguished signatures of reionization. Thereby, complementary to the constraints 
on the late stage of reionization 2; ~ 6 from observations of the most distant quasars absorption 
lines, etc, CMB provides a unique probe for the early stage of reionization. On the other hand, 
in order to interpret the observed spectra of CMB anisotropies and polarizations within the 
standard model, the reionization-induced modifications have to be taken into account properly. 
As is known, the reionization parameters could be entangled with the cosmological parameters, 
thus biasing our interpretation of CMB, and of reionization as well [9, 10, 11, 12, 13, 14]. In this 
regards, analytic studies can improve our understanding of CMB and reionization, even though 
the comparisons with the observed data need more accurate numerical calculations, such as 
cmbfast and CAMB [15, 16]. 

Two kinds of perturbations of the spacetime metric, i.e., density perturbation [17, 18, 19] 
and relic gravitational waves (RGWs) [18, 20, 21, 22], will effectively influence the CMB through 
the Sachs-Wolfe term [23] in the Boltzmann equation for photons. Although the contribution 
by density perturbation is dominant, RGWs give rise to a magnetic type of CMB polarizations, 
providing a distinguished channel to directly detect RGWs of very long wavelength [9, 11, 26, 27]. 
Moreover, RGWs have substantial contributions to large angular scales part of CMB spectra, 
where the impact of reionization is also dominant and cause bumps in the CMB polarization 
spectra for / < 10. Thus, in order to study reionization through the CMB, one has to take into 
account of the contribution of RGWs, or, vice versa. 

The analyses have been made towards CMB anisotropies and polarizations generated by 
RGWs [20, 21, 24, 25, 27, 28, 29, 30, 31]. In particular, by an approximate treatment of 
the time integration over the decoupling process during the recombination, Refs.[32, 33] have 
derived the analytic expressions of the CMB polarization spectra, C^^ and C^^. Recently, 
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extending the previous works, we have improved the time integration by a better approximation, 
and obtained the analytical expressions of all the four spectra, including Cf^ and Cj^ [34], 
which agree fairly with the numerical results up to a broader range of multipole moment / < 600. 
In that work the damping on RGWs due to neutrino free-streaming (NFS) has been included 
[35, 36, 37, 38], and its effects on the cross spectrum Cf^ have been demonstrated in details. 
In these analytical calculations, the reionization process has not bee included, which will be 
addressed in this paper. For the purpose of calculating the reionized CMB spectra C/^^, the 
reionization can be treated similarly to the decoupling, if the visibility functions for both processes 
are given. While the decoupling and its visibility function Vdirj) effectively distributed around 
z ~ 1100 have been better studied, the reionization is currently less understood, and is commonly 
modeled by its ionization fraction Xe{r]) as a function of time. We shall examine three possible 
simple reionization models with explicit Xe{rf), which, for a given value of the optical depth Kr, 
can be converted into its corresponding visibility function K(?7) effectively distributed around 
z ~ 11. The functions Vd{r]) and Vr{ri) are separately distributed, not overlapping, each of 
them can respectively be approximated by Gaussian type of functions, which are specified by 
their location, height, and width. In parallel, we will carry out, with approximation, the time 
integrations of Boltzmann's equation for the decoupling and reionization processes. The modes ai 
and (3i, respectively, for CMB temperature anisotropies and polarization, are obtained as analytical 
expressions. Each mode explicitly consists of two separated parts, one from the decoupling, and 
another from the reionization. Moreover, the optical depth Kr appears as the coefficients in ai 
and I3i in different combinations, and probabilistic interpretations are given. Besides reionization 
and decoupling, the result contains also other cosmological parameters for inflation that are 
contained in RGWs. Thus analytic studies on the reionization effects will be facilitated. 

In Section 2 we review briefly the result of RGWs spectrum h{u,r]) that will be used as 
the source for CMB anisotropies and polarization. In Section 3 three models of homogeneous 
reionization Xe{ri) are presented: one sudden and two extended. For each model the visibility 
function Vr{rj) and the optical depth function K,r{v) ^''s presented. In Section 4, by approximately 
carrying out the time integrations, the analytical expressions of ai and Pi are obtained. The 
resulting spectra Cf-^ are demonstrated. In Section 5, detailed analyses are made towards the 
reionization effects upon Cf--^ , three models are compared, and, in particular, examinations are 
made on the degeneracies of Kr with the normalization amplitude A and the initial spectral index 
Pinf of RGWs produced during inflation. The effect of reionization on the zero multipole analysis 
is addressed. The conclusion is given in Section 6. We use the unit in which c — h — ks — i ir\ 
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this paper. 



2. RGWs Spectrum 



The expansion of a spatially flat Universe can be described by the spatially flat {Qa + + 
Qr — 1) Robertson-Walker spacetime with a metric 

ds"^ = a^{r])[-drf + {5ij + hij)dx'dx% (1) 

where a{ri) is the scale factor, rj is the conformal time, and hij is the gravitational waves, taken 
to be traceless and transverse (TT gauge) = 0, and hij^j = 0. By the Fourier decomposition 

^..(^,x)=X:/^ef,.4-)(r?)e^'^- (2) 

for each mode k and each polarization a — (+, x), the wave equation takes the form 

/ife + 2-hk + k'^hk = 0, (3) 
a 

where the polarization index a has been skipped for simplicity, and the subindex k can be 
replace by k since the perturbations are assumed to be isotropic. The analytic solution of Eq.(3) 
has been given for the expanding universe with the consecutive stages: inflationary, reheating, 
radiation-dominant, matter-dominant, and accelerating, respectively in Refs.[38, 39, 40]. In our 
convention, 

ain) = am{r] - Vmf: V2<V< VE: (4) 

for the matter-dominant stage, and 

a{ri) ^ Inlrj - Val''^ , rjE<rj<'no, (5) 

for accelerating stage up to the present time r]o, where 7 ~ 1.044 for il^ = 0.75, and In = i/Hq, 
Hq is the Hubble constant. The normalization of a{ri) is chosen to be \r]o — ria\ = rja — r]o = I, 
where we have taken rjo = 3.11 to be the present time. Then, once the ratio Qa/O^ is specified, 
all the parameters will fixed: = ^if^G^^^^^^^ Ve ^ Va - Ce^ , Vm ^ Ve - ^Cb^ with 
(e = i^A/^m)^^^- The details have been explicitly demonstrated in our previous study of RGWs 
[38, 39]. 

When the NFS is included, a process occurred from a temperature T ~ 2 MeV during the 
radiation stage up to the beginning of the matter domination, the analytic solution hk{ri) has 
been given [34, 38]. The NFS causes a damping of the amplitude of RGWs by ~ 20% in the 
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frequency range (10~^^, 10^^°) Hz, leaving observable signatures on the second and third peaks 
of CMB anisotropies and polarization. So the RGWs damped by NFS will be used as a source 
in our calculation. As for other physical processes, such as the QCD transition and the 
annihilation in the radiation stage [37, 40, 41], they only cause minor modifications of RWGs on 
the small scales p > 10~^^ Hz, not being observable in the present large-scale CMB spectra, and 
will not be considered here. 

The solution hk{r)) depends on the initial condition during the inflation stage. We choose 
the initial spectrum of RGWs at the time rji of the horizon-crossing [34, 38, 39, 42] 

nu3/2 I. 

hiu,n,) = -^\hM\ = A{^r+f^-f, (6) 

where kn — 27r is the comoving wavenumber corresponding to the Hubble radius, A \s a k- 
independent constant to be normalized by the present observed CMB anisotropies in practice, 
and the spectral index Pin/ is a parameter depending on inflationary models. The special case 
of Pinf = —2 is the de Sitter expansion of inflation. If the inflationary expansion is driven by a 
scalar field, then the index (3inf is related to the so-called slow-roll parameters, r] and e [43], as 
Pinf = — 2 + (r] — 3e). Pinf is related to the spectral index ns of primordial scalar perturbations 
as ns = 2Pinf + 5. in literature, the RGWs spectrum is also written in the following form [1] [2] 
[44] 

Al(k)^Ar{lr = h\u,n,), (7) 

where the tensor spectrum index ht = '^{Pinf + 2) ~ without the running index, A-q is some 
pivot wavenumber, taken as fco = 0.002 Mpc~^ in our calculation, and the tensor spectrum 
amplitude At = 2.95 x 10^^ A(ko)r, where A{ko) is the scalar power spectrum amplitude that 
can be determined by the WMAP observations [1, 3, 4], and we take A{ko) ~ 0.8 accordingly. 
The tensor/scalar ratio r is model-dependent, and frequency-dependent [33, 45]. Recently, the 
5-year WMAP data improves the upper limit to r < 0.43 (95% CL) [8], and combined with BAO 
and SN gives r < 0.2 (95% CL) [5] [7]. In our treatment, for simplicity, r ~ 0.37 is only taken 
as a constant parameter for normalization of RGWs, except otherwise mentioned. 

The resulting functions hk{rj) and hk{rj) serve as the tensorial source to CMB anisotropies 
and polarization. Without reionization, only RGWs hk{rjd) and hk{rid) at the decoupling time 
rjd are relevant, contributing to the primary CMB spectra. When reionization comes, hk{r]r) 
and hk{r}r) at the reionization rjr contribute too, mainly contributing to the very large angular 
reionization bumps of CMB spectra. In Fig.l, hk{rid) and hk{rid), and hk{rjr) and hk{rir) are 
plotted. The right panel of Fig. 1 shows that, h{r]d) has the greatest amplitude around k ~ 25, 
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Figure 1: The RGWs hk{rid) and hk{rjd) at the decouphng and hk{rjr) and hk{rjr) at the 
reionization. 

forming a deep trough, whereas h{r]r) has the greatest amplitude around A; ~ 2, forming a deep 
trough. The left panel shows that both hk{r]d) and hk{rjr) have similar slope for small k. As we 
will see, these features of RGWs at r]d and at rjr are responsible for the profiles of CMB spectra 

3. Visibility Function 

in Basko and Ponarev's method, the Boltzmann equation of the photon gas for the k-mode 
is written as a set of two coupled differential equations [20, 21] 

ik + [ikji + q]^k = hk, (8) 

13k + [^k^x + q](3k = qGk. (9) 
where (5). is the linear polarization contributed only by linearly polarized CMB photons, 

ttfe = Cfe - (10) 

is the anisotropy of radiation intensity contributed by both unpolarized (natural light) and polar- 
ized CMB photons, ji — cos 9, q is the differential optical depth, and 

Gkiv) = ^ £ v[(i + ^^yPk - ^(1 - (11) 

In the following, we omit the subscript k for simplicity of notation. The formal solutions of 
Eqs.(8) and (9) at any time rj can be written as the following time integrations [33, 34]: 

e(77)= r h{r]')e-^^''''''^e''"'^'^'-'^Urj', (12) 
Jo 

p(r))^ r G(r)')q(r)')e-''^'^''^">e'''^'^'^'-'^^dr)', (13) 
Jo 

where 

rv' 

K{r]',r])= / qdr] = k{7]) — K{rj') (14) 

Jt) 

with the optical depth given by 

K{ri) = K{rio,rj) = q{rj')dri' (15) 

Jri 
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from the present time 7]o back to an earlier time 7], such that 

= (16) 

The CMB anisotropics and polarization are usually expressed in terms of their Legendre compo- 
nents 

m = \f_Ji^i{v,iAPi{iA, (17) 

Pi{v) = lj'^dpi/3{rj,pi)Pi{pi), (18) 
where Pi is the Legendre function. By the expansion formula 

oo 

e-/^ = ^(2/ + iyj,(x)PK//) (19) 

and the ortho-normal relation for the Legendre functions, the components at the present time 
770 are given by the following 

6(%) = r e-''^''^Hv)ji{k{v - Vo))dV: (20) 

PiiVo) = r G{ri)V{ri) ji{k{rj - rio))drj, (21) 
Jo 

where 

V{r]) = q{r])e~<'^^ (22) 

is the visibility function. As one sees, to analytically carry out the integrations in Eqs.(20) and 
(21), one needs the explicit expression of e~'^'^'^'> and V{vi), which are determined by the whole 
history of ionization. In the following we will give approximate formula of both functions. 

V{r]) has the meaning of the probability that a CMB photon reaching us today was last 
scattered by free electrons at the time 77. Without the reionization, V{r]) would have only one 
sharp peak around z ~ 1100 for the decoupling, and satisfies the normalization condition 



rm 

/ Vir])dr] = 1. (23) 
Jo 



When the reionization is included, V{ri) will have, around z ~ 11, another peak. If the universe 
was reionized twice, say at 2; ~ 6 and z ~ 16, [46, 47, 48], V{ri) would have double peaks 
for reionization. We consider only the case of a single reionization in this paper. Then, as a 
function of t], V{r]) is mainly distributed around decoupling and reionization, and is effectively 
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Figure 2: The visibility function Vd(77) for the decouphng around z ~ 1100. Both the 
analytic and the fitting by two half-gaussian functions are shown. 



vanishing in the region far away from the peaks, as shown in the Panel (d) in Fig. 6. Thus the 
time integration of Eq.(23) can be practically split into two parts 



JO 



VMdv+ Vr{v)dv = l, (24) 



'0 J Vsplit 

where Vd{r]) and K(^) are the portions of V{rj) for decoupling and reionization, respectively, 
and rispiit is some point between decoupling and reionization with V{r]spiit) — 0. In calculation 
we can take, say, tjspiu — 0.297 corresponding to a redshift z ~ 100. In Eq.(24), Jq'^''^ Vd{r])dr] 
is the area covered under the curve of Vd{7]), and stands for the probability that a photon was 
last scattered during the decoupling. Similarly, J^°^^^^ Vr{r])dr] is the probability that a photon 
was last re-scattered during the reionization, i.e., the amount of CMB photons out of the total 
that are rescattered. According to Eq.(24), their sum is constrained to be unity. This has a 
physical interpretation: more CMB photons are last scattered around ~ r]r, less will be last 
scattered around ~ rid- During reionization the intrinsic anisotropies of this portion of CMB 
photons were washed out, and new polarizations were generated on large angular scales. As we 
will see, J^°^^.^ Vr{ri)dr] depends essentially on the optical depth up to the reionization. 

Now let us specify the visibility functions Vd{r]) and K('?)- First the decoupling process is 
better understood, whose Vd{r]) has been given explicitly, which depends the baryon fraction VLb 
[19, 49, 50]. As a function of time, the profile of Vd{rf) itself looks like a sharp peak around the 
decoupling z ~ 1100. Thus, when it appears as a factor of the integrand in the time integration 
(21) for the polarization A(^o). 't actually plays a filtering role: only the narrow time range 
around the decoupling contributes substantially to the integral of Eq.(21). To facilitate analytic 
calculations of CMB polarization, Vd{ri) has been approximated by the following two pieces of 
half gaussian function [33, 34] 



Vd{^) 



n^<^)exp(-^K (^<%), 



(25) 



where rid is the decoupling time, which is taken rid — 0.0707 corresponding to a redshift Zd — 
1100, A77di = 0.00639, Ajy^s = 0.0117, and {Aridi + Arid2)/2 = Arid = 0.00905 is the thickness 
of the decoupling. Eq.(25) improves a single gaussian function [32] by ~ 10% in accuracy and 
at the same time allows an analytic treatment of the CMB polarization spectrum. We have 
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Figure 3: The three models of reionization with a fixed optical depth Kr = 0.084. For each 
Xe{ri) given in Eqs.(26), (27), and (28), the functions qr{i])j i^rij])i ^^^d K-(^) are calculated 
according to the formulae in Eqs.(29), (30), and (31), respectively. 



checked that the errors between Eq.(25) and the numerical formulae given in [19, 50] is very 
small, < 3.9% in the whole range. The coefficient V{rid), as the height of Vd{'r]), also depends 
on the reionization through the normalization in Eq.(24). The analytic Vdijf) with VLb = 0.046 
and its fitting are shown in Fig. 2. 

Next, understanding of the reionization as a physical process is still underway, and various 
tentative models have been proposed for it. Spatially, the reionization might have occurred 
inhomogeneously [51, 52, 53, 54, 55], resulting in modifications on the small angular scales 
part of CMB spectra. Models of double reionization [46, 47], or its variants, such as peak-like 
reionization [56], have also been proposed. In the following, we will work with three simple 
homogeneous models, whose ionization fraction Xe{rj) are explicitly given. One is the sudden 
reionization model with 

where r]r is the reionization time. For concreteness of illustration, in our calculation we take 
r]r = 0.915, corresponding to the redshift Zr = 11. This is the simplest model often used 
in the literature. But there are accumulating evidence that the reionization is an extended 
process, stretching from 2; ~ 6 up to 2; ~ 11, even up to as early as 2; ~ 20 [8, 57, 58]. For 
instance, studies of Lya Gunn-Peterson absorption [59] indicate a rapid increase in the ionized 
fraction of the intergalactic medium at a redshift lower than Zr — 6. On the other hand, the 
WMAP observations of CMB found a much earlier reionization, 2;^ = 17 ± 5 by WMAP 1-yr 
[1], Zr = 10.9l|^ by WMAP 3-yr [3], Zr = 11.0 ± 1.4 (68% CL) by WMAP 5-yr [8], and 
Zr = 10.8 ± 1.4 by WMAP 5-yr combined with SN and BAO [7] [5]. One extended reionization 
model is the 77-linear reionization with 



0, for r] < rjri 
for 77,1 < 77 < Vr2, (27) 

1, for r] > r]r2- 



where rjri and 77^2 are the beginning and end of reionization. For instance, one can take rfri = 
0.685 and 77^2 = 1.20712, corresponding to Zri = 20 and Zr2 = 6, respectively. This model is 
closer to the result of WMAP 5-yr fitted by the two step reionization [8]. Another extended 
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reionization model is the 2;-linear model with [55]: 



0, 
1 - 



Z — Zr2 



for Z > Zrl 
for Zrl > Z > Zr2, 



Zrl-Zr2 ' 

1, for Z < Zr2- 



(28) 



For Zrl — 20 and Zr2 — 6, one has Xe{z) = 1 — — 6)/14. The ionization fraction Xe{r]) for 
these three reionization models are comparatively shown in Fig. 3. 

Given Xe{r]) in the above three models, the differential optical depth for reionization can be 
directly calculated by the formula [19, 55, 60]: 



(29) 



where the constant = (1 - ^)^h£c£T^ y ~ 0.23 is the primordial helium fraction, ar 
is the cross section of Thompson scattering, nip is the mass of a proton. For Jl^ = 0.045, 
Cc — 0.142 X 10~^^ m~^. Since the value of Yp from observations has considerable large error 
bars [61], in our treatment Cc is allowed to vary slightly around this value. From Eq.(15) follows 
the optical depth for reionization as an integration 



and, from Eq.(22) follows the visibility function for the reionization. 
For instance, for the sudden reionization model, one easily obtains 



(30) 



(31) 



«r(?7) = (V-Vm) -(VE-VmY 



3 



-I 



H 



(32) 



27 + 1 

where all the parameters have been given bellow Eq.(5). For a reionization model, the most 
important quantity k,,. = Kr{r]b) is the value of the optical depth from r]o back up to some 
time rjb before the reionization, where QriVb) is practically vanishing. For example, one can take 
Vb = Vspiit- In practice, one can conveniently take r]i, — r]r for the sudden model, and take 
77b = 0.5 for the ?7-linear and 2;-linear models. Kr is an integral constraint on the reionization 
history. On the observational side, based upon treatments of a sudden model, WMAP 1-yr 
gives Kr = 0.17 ± 0.04 [1], and WMAP 3-yr gives Kr = 0.09 ± 0.03 [3], and WMAP 5-yr gives 
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Figure 4: The r/-linear reionization model with Kr = 0.084. The sohd hnes are the calculated 
results. The dashed lines are the fitting by two half Gaussian functions in Eq.(33). 



Figure 5: The ^-linear reionization model and its fitting. 



Kr = 0.087 ±0.017 [8], and WMAP 5-yr combined with SN and BAO yields Hr = 0.084 ±0.016 
[5, 7]. To be specific in calculation, we will take the value = 0.084 for all three reionization 
models in this paper, except when it is mentioned otherwise. However, note that, for extended 
reionization models, one should be careful in applying the WMAP observed value of Kr, as it 
is obtained by using a sudden model. For the r^-linear model with Xe{r]) given in Eq.(27), one 
uses the formulae of Eqs.(29) (30) (31) to compute qr{r]), K,r{r}), Vr{ri). For the z-linear model 
with Xf.{r]) in Eq.(28), one does similar computations. The resulting qr{r]), Kr{r]), and K(^) foi" 
these three models are plotted in Fig. 3. 

The value of optical depth Kr determines the area S^°^^.^ Vr{r])dr] introduced in Eq.(24). For 
a fixed Kr = 0.084, the integration of Eq.(31) yields J^°^^.^ Vr{rj)drj = 0.0795 in the sudden 
model, Vr{r])dr] = 0.07953 in the 77-linear model, and 1^1°^^.^ Vr{r])dr] = 0.07973 in the 
2;-linear model, respectively. So two gradual models have slightly larger area than the sudden 
model. Besides, our computations also show that a larger Kr yields a larger f^°^^.^ Vr{r])dr] and a 
smaller Jj'^*'*'* Vd{r])dr] due to Eq.(24), meaning that a CMB photon reaching us was more likely 
last scattered at reionization. As we shall see explicitly, for CMB spectra, this will enhance the 
reionization bumps on large scales and reduce the primary peaks due decoupling. 

To facilitate analytical calculations of CMB polarization, similar to the treatments of Vd{ri) 
for the decoupling, Vr{ri) can be also approximated by some fitting formula. For the 77-linear 
model, it is fitted by the following two pieces of half Gaussian functions 



where Ar/^i = 0.147, Ar]r2 = 0.425, Ar]r = {Ar]ri + Ar]r2)/2 = 0.286, and r]r = 0.935 
{zr = 10.5). It is plotted in Panel (c) of Fig. 4 under the requirement that it gives the same 



Eq.(33) but with the parameters Ar]ri = 0.100, Ar]r2 = 0.366, Ar]r = {Ar]ri + Ar]r2)/2 = 0.233, 
and 7]r = 0.855 {zr = 13). It is plotted in Panel (c) of Fig. 5. Here for the two extended models, 
the value of rjr has been taken to correspond to the maximum of K(^)- For the sudden model. 




iV < Vr), 
iV > Vr), 



(33) 



area /^^"^^.^ Vr{ri)dri as the calculated one. For the 2;-linear model, the fitting formula is similar to 
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Figure 6: The sudden reionization model and its fitting. Panel (c) shows that the fitting 
Vr{i]) by Eq.(34) has large errors to the calculated one. The evolution history of V{rj), 
including both reionization and decoupling, is sketched in Panel (d). 



it can be fitted by a half piece of Gaussian function 

{0, (for 1] < r]r), 

nry.)exp(-i^), {iorv>Vr), ^^^^ 

with the width Ar]dr — 0.247, plotted in Panel (c) of Fig.6. The half-gaussian fitting of K(^) fof 
the sudden model is not as accurate as those for the two extended models. It should be expected 
that in the sudden model the analytical CMB spectra C/*-^ based on its fitting formula (34) is 
not as good as those in the two extended models. 

We mention that, given a fixed Kr, the respective height V{r]r) in Eqs.(26), (34), and (33) 
are also determined automatically. From these fitting Vr{r]), one can convert it to obtain the 
corresponding optical functions 

e->^r{v) = 1 - / Vr{r])dr], (35) 

Jt) 

Kriv) = - In (l - Vr{v)dv^ , (36) 

qriv) = 7 ^ V • (37) 

It should be mentioned that the approximate fitting of K(^) by Eq.(33) underestimates the 
value of Vr in the range 77 > 77^ by ~ 9.1%. For the linear model, the fitting by half Gaussian 
functions underestimates the value of Vr in the range 77 > 77^ by ~ 8.6%. However, this kind 
of error of the fitting can partially compensated in treating the damping factors occurring in the 
time integration of the polarization mode, as will be given in the following. The gaussian fitting 
of Eq.(34) for the sudden model is included only for illustration purpose, as its error is larger 
than the two extended models. 

4. Spectra of CMB Anisotropies and Polarization 

By applying the same kind of approximate integration technique as in Refs.[33, 34], up to 
the second order of a small l/q^ in the tight coupling limit, the function G{r]) in Eq.(ll) can be 
written as 

^(77) r hi-q') e-To^^'^'^-To^^Uri', (38) 

10 Jo 
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and the integration of polarization mode in Eq.(21) is written as 

PiiVo) = -^t' r dvV{r))h{r))ji{kir) - rjo)) P dr)'e-^o-(r,')-{-o-(v) _ (39) 
ID Jo Jo 

Since the visibility function ^(77) for the whole history consists of two effectively non-overlapping 
functions, Vd{r]) and Vr{r]), the 77-time integration Jq° drj in the above is naturally split into a 
sum of two integrations: 



1 fVsplit . fV 

PliVo) = -j^i^ dr]Vd{v)h{ri)ji{k{r] - rjo)) dr]'. 



g 10 "-V t 10^ 



i dr)Vr{r))h{r])ji{k{r]-r]o)) l dr)'e-^''^'^'^-^''^'^K (40) 



Vsplit 



One defines the integration variable x = K,{ri') / K,{rj) to replace the variable 77' in the above. Since 
Vd{r]) is peaked around rjd with a width Arjd, and, similarly, Vr{r]) is peaked around rjr with a 
width Arjr, one can take drj' ~ —Arjd^ and drj' ~ —Arjr^ as approximation, respectively. 

A(^o) = -TTT^'Ar^, / drjVd{v)h{rj)ji{k{rj - Vo)) / -e-To«W-To«W 
lU Jo Ji a; 

-^^Afjr r dr]Vr{v)h{r])mr] - r]o)) H -e-^o<v)^-^o<^\ (41) 

iU JVsplit Ji ^ 

For each term in the above, the r^-time integration can be dealt with, using the same kind of 
treatment as in Ref.[33, 34]. For the decoupling one has 

split . . rVsplit 

/ dr)Vd{v)Kv)jiiKv - Vo)) ^ Dd{k)hivd)ji{kivd - Vo)) / drjVd{v), (42) 
Jo Jo 

where the damping factor for the decoupling is given by the following fitting formula 

Dd{k) = —[e-'^kAr,,,]^ ^ g-cCfcAr^.^)"]^ (43) 

which can be simplified by 

Dd{k) = 1.46-"^'=^'''')', (44) 

with c and b being two fitting parameters. For CMB spectra without reionization, it has been 
shown in Ref.[34] that both damping factors in Eqs.(43) and (44) c ~ 0.6 and b ~ 0.85 give a 
good match with the numerical result by CAMB [16] over an extended range I < 600, covering 
the first three primary peaks, and the error is only ~ 3%. 
Similarly, the 77-time integration for the reionization is 

rvo . . rvo 

/ dr]Vr{r])h{r])ji{k{r} - r)o)) ^ Dr{k)h{r]r)ji{k{r)r - r)o)) dr)Vd{r]), (45) 

J^split J^split 
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where the damping factor for the extended models is taken to be 



2 



or for the sudden reionization 



Dr{k) 



(46) 



(47) 



Here the parameter c and b in Eqs.(46) and (47) for reionization could take values different from 
those for decoupling. For simplicity, we let them take the values that are the same as in Dd{k). 
Guided by the error estimation for the decoupling case, we can only estimate the errors due 
to Dj.{k) in Eq.(46) for the two extended models upon the reionization bumps of polarization 
spectra to be < 10%, the same order of magnitude as those of the fitting K-(^) in Eq.(33). 
Substituting Eqs.(42) and (45) into Eq.(41), and performing the integrations / drj first, 



rVsplit 3 / N 7 / ^ 3 

/ d77Vd(77)e-^o«('?)^-To«('?) = / dKe-To 

Jo J Kr 



1 n ir\^ 



10 ^ 10 



rvo __3_ 
/ dr]Vr{r])e " 

•'Vsplit 



due lo'* = 



1 



10 ^ 10 



1 — e (ts+^^)'*'' 



(48) 
(49) 



one finally obtains the expression of the polarization mode as a sum of two parts 

PiiVo) = [Ai{Kr)Dd{k)Ar]dh{r]d)ji{k{r]d - rjo)) 

+A2{Kr)Dr{k)Ar]rh{r]r)ji{k{r]r - rjo)) 

where the re^-dependence coefficients 

/" 

Ai{Kr) = 



dx 



X 



A2{Kr) 



17 
10 

dx 



"( 10 + in^)'*'' 



10 



X] 



^(lO + 10^) 



(50) 



(51) 



(52) 



both being independent of the wavenumber k, and the sum is Ai(/tr)+A2(/tr) = In ^ ~ 1.116 



independent of Kr- If one sets A2 = and Ai = ^In y, Eq.(50) reduces to exactly that of the 
non-reionization case [33, 34]. Actually, after the sum is normalized to unity, the two coefficients 
have the physical meaning: 

Al{Kr) 



ai(Kr) 



10 In 20 

17 3 



(53) 



is the probability that a polarized photon we perceive was last scattered during the decoupling 
epoch, and 



_ A2{Kr) 
a2[Kr) — 10 , 20 
17 3 



(54) 
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Figure 7: The normalized coefficients ai^Hr) and a2(Kr) of the polarization Pi. A larger 
Kr yields lower ai{Kr) and higher a2{K,r), i-e., will yield lower primary peaks and higher 
reionization bump in Cp^ and C^^ . Also plotted are the coefficients e^"'' and (1 — e"'*'') of 
the temperature anisotropies Notice that ai(/tr) and 02 (/tr) vary with Kr more drastically 
than e"'^'' and (1 — e"'*''), respectively. 

is the probability that a polarized photon we perceive was last scattered during the time interval 
from the beginning of reionization up to the present time 770- It is found that ai{K,r) is a decreasing 
function of Kr and a2{Kr) is an increasing one, as shown in Fig. 7. Therefore, if more CMB 
photons are scattered by the free electrons during the reionization, the optical depth Hr acquires 
a larger value, giving rise to a higher coefficient ^2(«^r) and, at the same time, a lower coefficient 
Ai{Kr). The Ai{Kr) part in Pi from the decoupling will give rise to the primary peaks of Cf^ 
and C[^^ , and will be prominent on small angular scales with / > 100. The A2{Kr) part from 
the reionization will be dominant on large angular scales and will yield the reionization bumps of 
Cf ^ and CP^ around / < 10. 

The analytical expression (50) has the merit that effects of relevant physical elements upon 
the polarization have been explicitly isolated and displayed. The K^-dependence of (3i is attributed 
to the coefficients Ai{Kr) and A2{Kr), which determine the relative heights of the primary peaks 
and the reionization bump. Other effects of reionization is encoded in the factor Dr{k)Arir- The 
effects of decoupling are absorbed in Dd{k)Arid. The effect of RGWs upon the polarization are 
given by the time derivatives h{ria) at rj^ and h{r]r) at rjr, which not only contain the cosmological 
information, such as inflation and NFS, etc., more importantly, but also determine the overall 
profiles of Cf^^ and C^^, such as the locations of peaks and troughs, and of bumps. The factors 
ji{k{r]d — rjo)) and ji{k{r}r — rjo)) just play the role of conversion from the wavenumber /c-space 
into the multipole Z-space. 

To calculate the temperature anisotropies, we need to evaluate ^i in Eq.(20), which contains 
the factor e~'^^'^\ This also needs to be dealt with properly. As shown in Fig.6, the factor e"**^^^ 
has two steps, one at the decoupling 77 = rj^, and another at 77 ~ 77^ caused by the reionization. 
It can be approximated by the following two-step function 



iv < Vd); 

e""-" {Vd<V<Vr); (55) 

1 {r]r<r]< Vo), 



and its reionization-relevant part e '^'^'^^'^ is between (77^,770). By Eq.(35), e '^'''^^^ is the integration 
of Vr{rj) from 77 to 770, determined by the area under the curve of K(^). riot very sensitive to the 
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detailed shape of Therefore, the approximate formula (55) will be used for the three models 

of reionization, with their respective values of r]r. Note that Eq.(55) tends to overestimate the 
contribution of the reionization to the integration, since e"'^*^'') shown in Fig. 6 increases gradually 
from e"'*'" at r]r up to 1 for 77 ^ rjr, instead of instantaneously jumping up to 1 at rjr- To 
compensate this overestimation, in actually calculating Ci{rjo) in the linear model, we may use 
the value of slightly greater than 0.935. But this adjustment of the time T]r does not apply 
to (3i{r)o) in Eq.(50). Substituting Eq.(55) into Eq.(20), the integration for is split into two 
terms 

^i{r]o) ~ i' / e-'''^h{rj)ji{k{rio - rj))drj + / h{ri)ji{k{rjo - rj))drj. (56) 

Following the similar treatments in [45, 34], each term is integrated by parts, yielding the following 
approximate expression 

ii{r]o) = -i' [e''''''h{r]d)ji{k{r)Q - Vd)) + (1 - e-^^)h{r]r)jiik{r)o - Vr))] , (57) 

where the first term is generated by h{rid) at the recombination and the second term is due 
to h{r]r) at the reionization. Eq.(lO) then yields the mode of CMB temperature anisotropics 
<^i{Vo) — 6(^0) — Piivo)- In fact, ai{rjo) is essentially contributed by ^^(770) since the amplitude 
of 6(^70) is about two orders higher than that of A(^o)- Writing down explicitly, one has the 
approximate, analytic expression of the mode of CMB temperature anisotropics, including the 
reionization. 



-i^jiiHvo - Vr)) 



e "^hirjd) - ^Ai{Kr)Dd{k)Ar)dh{rjd) 
(1 - e-^'-)h{r]r) - ^A2{Kr)Dr{k)Ar]drh{r]r) 



(58) 



In this expression, the first term containing h{r]d) and h{7]d) is brought by the decoupling and 
responsible for the primary peaks, whereas the last term containing h{r]r) and h{rir) is brought 
in by reionization and prominent on large angle scales with I < 10. When one sets Ai = 1, 
A2 = 0, and e^'^'' = 1, Eq.(58) reduces to the results for the case without reionization [34]. The 
K^-dependence of ai is mainly attributed to the factors e""'" and (1 — e^**' ), while the portion 
containing Ai^Kr) and A2{K,r) is the subdominant (3i. By Eq.(35) and the definition of Kr, on 
has 

Vriv)dv, (59) 



Vb 



which has a physical interpretation: the probability of a CMB photon being last scattered during 
the earlier epoch before the reionization. Since e"'*'^ < 1 for Kr > 0, it will cause a slight decrease 
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in the amplitude of the temperature anisotropies, as demonstrated in Eq.(58). Correspondingly, 
the factor (1 — e"'*'') in front of h{r)r) is 



rvo 

= / Vr{r])dr], 

Jrib 



(60) 



recognized as the probability of a CMB photon being last scattered during the time interval from 
the reionization up to the present time 770- These foregoing probabilistic interpretations have the 
parallels in the case of CMB anisotropies generated by scalar perturbations, where reionization 
also brings about a similar exponential factors e~'^'^ in the temperature anisotropies, and a physical 
illustration on its appearance is given in Ref.[62]. It should be mentioned that the probabilities in 
Eqs.(59) and (60) are respectively different from the normalized a-i{Kr) and a2{Kr), the latter are 
for the polarized photons. Moreover, as shonw in Fig. 7, ai{Kr) decreases with Kr much faster 
than e"'^'^ does, and a2(Kr) increases much faster than (1 — e"*^"^). In this sense, the polarization 
l3i{rio) is more sensitive to Kr than the temperature anisotropies Q!;(?7o). Therefore, one may say 
that the polarization spectra Cf^ are Cf^ are more sensitive probes into the reionization than 
the temperature anisotropies spectrum Cf^. 

With ai and fii being ready, one can compute straightforwardly the CMB spectra caused by 
RGWs. The detailed derivations have been demonstrated in Refs. [27, 33, 34]. In particular, 
some minor misprints of the coefficients in Ref.[27] have been pointed out and corrected in Refs. 
[33, 34]. The temperature anisotropies 
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(62) 



the magnetic type of polarization 
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Substituting ^^(^^o) and Pi{rio) into Eqs. (61), (62), (63) and (64) yields the analytical expressions 
of the spectra of CMB with the modifications of reionization: 



TT 
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+PTl{k{V0 - Vr)) 
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(1 - e ^^-Mrjr) - —A2Dr{k)Ar]rh{rir) 
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In the above integrations , the projection factors are defined as: 

jl-2{x) , 2ji{x) , Jl+2{x) 



Pti{x) = 
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(70) 



(71) 



We apply these formulae to the three reionization models, respectively, and plot the spectra 
C^-^ . The reionized spectra C^-^ are plotted in Fig. 8 for the three models of reionization, in 
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Figure 8: The spectra in three reionization models. The parameters c = 0.65 and 

b = 0.85 in Dj.{k) are taken. The optical depth k,,. = 0.084 and the ratio r = 0.37 are taken. 
The numerical result is obtained with the same set of parameters, using CAMB [16]. 

which we also plot the numerical spectra from the CAMB Online Tool for a comparison [16]. 
Both the analytic and numerical computation use the same set of parameters Kr = 0.084 and 
r = 0.37. On large scales / < 600 our analytical Cf^ and C/^^ agree with the numerical ones. 
For the two extended models, the error is ~ 3% for the primary peaks, and the error is estimated 
to be < 15% for the reionization bumps as superposed from that of decoupling ~ 3% and that 
of reionization ~ 10%. Notice also that the analytical C^^ and C^^ in the sudden model have 
reionization bumps too low. This has been expected, since the half-gaussian fitting formula (34) 
is poor. The analytical Cf^ and Cf^ are close to the numerical ones on smaller scales I > 20, 
but have obvious departure from the numerical ones on very large scales I < 10. This implies 
that the approximation of temperature anisotropics in Eq.(57) is poor for small multipoles 
I < 10. In the following we focus only on the two extended models and examine the impact of 
reionization through the analytical spectra C/^^. 



1. The most prominent modification due to the reionization is that it enhances the low-/ parts 
of the spectra, forming a reionization bump at / ~ 5 for C^^ and C^^ , respectively. The position 
of this bump is a reflection of the horizon scale at reionization, whose corresponding angular scale, 
/ ~ 5, is much larger than / ~ 100 of the primary peaks at the photon decoupling. As pointed out 
earlier, the profiles of C/^^ are determined by the profiles of RGWs at the decoupling and at the 
reionization as well. In particular, the reionization bumps are generated by h{rjr), and the primary 
peaks and troughs are due to h{rjd) and h{r]d). This correspondence is clearly demonstrated by 
Fig. 9, in which the left panel plots C^^ and C^^, as well as h{r]d) and h{r]r) in one graph, and 
the right panel plots C^^, as well as h{r]d) and h{r]r) in one graph. This correspondence can be 
further explained by the following analysis. The respective projection factors, Pti, Peu and Pbi 
as the integrands of C/^^ are made up of the spherical Bessel's functions, ji{x), which is sharply 
peaked around x ~ Z. Subsequently the projection factors as functions of k are sharply peaked 
around 
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Figure 9: The correspondence of the profiles of Cf'^ with that of RGWs at the decoupling 
and at the reionization. h{ri,) is responsible for the bumps of Cf^ and C^^ around / ~ 5, 
while h{r]ii) is responsible the primary peaks and troughs for / > 100. For Cf^, ^(%) and 
h{r]r) have a similar slope at I < 10, and their superposition does not form a prominent 
bump of Cf^. 

respectively. Consequently, the spectra as integrations over k will receive main contributions 
from the integration range k ~ Z/770 to the primary peaks and from k ~ l/ivo — Vr) to the bump, 
respectively [33]: 



C^: Cr oc AlDlik) h{ria) " , + AiD'^{k) h{rir) [ . (74) 



2 



oc e-^- |M%)lL/.o + (1 - ^-^-y \hiVr)\Lin„.) , (75) 

C^^ oc Aie-^^Da(k)h(r]d)h(r]d)kr.i/r,o + ^2(1 - e-''^)Dr{k)h{r]r)h{Vr)k^i/{vo-vr)- (76) 

According to Eq.(74), the locations of the primary peaks Cf^ and Cj^^ are mainly determined 
by the \h{rid)\'^ -term, and those of the reionization bumps are determined by the \h{rir)\'^ -t&rm. 
However, the spectrum Cf^ does not have a prominent bump around / ~ 5. This is because both 
|^(%)P and l^(^r)P have a similar slope around there, and their superposition only enhances 
the spectral amplitude, not forming a bump. These are illustrated in Fig. 9. 

2. The reionization bumps in the polarization spectra depend on the detailed reionization 
history. Cf'^ and Cf^ for a fixed value of the optical depth Kr — 0.084 in the two extended 
models are shown in Fig. 10. The bumps in the 77-linear model are located at a slightly larger 
angular scale (smaller I) than that in the 2;-linear model. This is because we have assigned a 
greater rjr — 0.935 in the 77-linear model than that rjr — 0.855 in the ^-linear model, so its bump 
is located at a slightly smaller I ~ k{riQ — rjr). Notice also that the 77-linear model produces 
higher bumps than the 2;-linear model. This is due the fact that the 77-linear model has a greater 
width Ar]r — 0.286 than that A77r = 0.855 in the 2;-linear model. Thus we conclude that the 
location of bump is quite sensitive to the the reionization time 77^, and the height of bump is 
sensitive to the width A77J. of reionization process. This feature is helpful to probe rjr and Arjr 
only if observational data on the bumps are accurate enough. However, when we let the two 
models to have the same set of parameters r]r and Ar]r, their reionization bumps predicted by 
our analytical formulation are very similar. The lesson is that the bump is an integrating result 
from the ionization fraction Xe{r]), and, in this regards, two different reionization histories via 
Xe{ri) can lead to similar bumps, as long as they have similar Vr{ri) [63, 64]. 
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Figure 10: Cf^ and C^^ in the extended reionization models. The two models yield different 
reionization bumps at Z ~ 5 since they are assigned with different values of rjr and Arjr. 

3. The overall profiles of CMB spectra are very sensitive to the optical depth of reionization. 
In particular, Kr is strongly degenerate with the normalization of the amplitude of primordial 
fluctuations, and this fact has been one of main difficulties to probe the details of reionization 
process [9, 10, 11, 12, 13, 65, 66, 67]. It should be emphasized that the reionization does not 
change the primordial amplitude A of RGWs in Eq.(6), which is implicitly contained in h{r]) and 
h{r]). The impact of Kr is through the coefficients Ai{K,r) and A2{K,r) in A Eq.((50)), as 
well as the coefficients e"'*'^ and (1 — e~'^'^) in ai in Eq.(58). The main features of the Kr — A 
degeneracy are clearly revealed by the analytical estimations in Eqs.(74), (75), and (76). 

For instance, look at Eq.(74) for C^^ and Cf^. A larger Hr gives smaller Ai and larger A2, 
leading to lower primary peaks and higher bumps of Cf^ and C^^, as illustrated in Fig. 11. But, 
this lowering of primary peaks can be compensated by a choice of a higher amplitude normalization 
A, which enhances the amplitude of h{r]d), resulting in the unchanged term Al{Kr)\h{r]d)\'^ , so 
that the primary peaks remain the same. This is the Hr — A degeneracy. Similar degeneracy in 
Cf"^ and Cf^ are also understood by Eq.(75) and Eq.(76). 

The Kr — A degeneracy can be broken. Again, take C^^ and C^^ as an example. While 
a larger k,. and a higher A can yield the unchange primary peaks, the reionization bumps get 
doubly enhanced, since the bump term Al{Kr)\h{r]r)\'^ in Eq.(74) gets doubly enhanced. This 
suggests a possible way to break the degeneracy. Eq.(74) tells that the relative height of the 
primary peaks and the bump is given by 



primary peak amplitude 

oc 

bump amplitude 



Al^Kr) 
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2 

h{r]r) 



(77) 



For any given RGWs, the ratio h{7]d) / h{Vr) is independent of A and completely determined, 

so one has ^ 

primary peak amplitude / Ai(/t,.)\ 

bump amplitude \A2{Kr) J 

This ratio only depends on the value of Kr and is not sensitive to the details of a reionization 
model. Therefore, using this ratio of heights, one can infer the value of Kr from the observational 
data of Cf^ and Cf^, thus breaking the degeneracy. 

4. The primordial fluctuation spectral index introduced in Eq.(6) is a very important 
parameter for inflationary models. Given a normalization A of the RGWs amplitude, a large 
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Figure 11: The Kr — A degeneracy. A larger value of k,,. enhances the bumps at / ~ 5 and, 
at the same time, reduces the primary peaks of C^^ and C^^ . The plot is made for the 
^-linear model. This degeneracy behavior also exists in the ?7-linear model. 

Figure 12: The Kr — Pinf degeneracy. Although the bumps and the 1** primary peaks are 
degenerate, the 2"*^ and S^'' primary peaks show clear departure. The plot is made for the 
2;-linear model. 

Pinf tilts the spectrum h{u,r]i), in such a way that RWGs is more strongly enhanced on smaller 
scales [39, 38]. The RGWs-generated spectra C^^ are subsequently tilted in the same way 
[33, 34]. Therefore, a larger P^f brings about a similar effect on Cf^ as a smaller Hr does, 
leading to certain bias in determining Hr [67, 68, 69, 70, 71]. Take Cf^^ and Cf^ as an example, 
for which the effect is more prominent. Fig. 12 shows that, for the 2;— linear model, the case 
{Pinf — — 2.02,«;r = 0.106) and the case {Pinf — — 2.10, = 0.084) yield almost overlapping 
curves of the bump and the 1** primary peak as well. 

The Kr — Pinf degeneracy can also be understood by the analytical estimation in Eq.(74). 

2 

While a large Pinf enhances h{'rjd) on small scales, a large Kr suppresses Ai{Kr), resulting 

2 

in an unchanged combination Ai(fi;^)^ h{Vd) ^or the primary peaks. But this degeneracy is 
clearly broken from the 2"*^ primary peak on. This is because the K^-induced change in Ai{Kr) is 
scale-independent, whereas the /5i„/-induced change in \h(r]d)\ depends on the scale. Therefore, 
one expects that data of the smaller scale C^^ and C^^ will be helpful in breaking the Kj. — Pinf 
degeneracy. Note also that the principal component method developed in Ref.[71] can protect 
the bias of Kr caused by Pinf. 

5. Although the magnetic type of polarization Cj^^ is thought to be a "smoking gun" of 
detection of RGWs, its detection is not done yet, which may be accomplished by a future CMBpol 
experiment [72]. The 5-year WMAP [5, 6] has given the observed cross-spectrum Cf^, which is 
negative (anti-correlation) in a range / ~ (50,220). Yet this observed Cf^ is a superposition of 
contributions by both scalar perturbations and RGWs. In order to extract the traces of RGWs out 
of Cf^ , one still needs to disentangle the contribution by RGWs from the total. In the so-called 
zero-multipole method [25, 45, 73, 74], one examines the impact of the tensor/scalar ratio r 
upon the zero multipole Iq around ~ 50, where Cf^ first crosses the value and turns negative. 
However, there are other factors that can influence the value of Iq. The variation of Iq caused 
by NFS has been estimated to be small AZ < 4 [34]. Here the reionization is another important 
factor that brings about a change of Iq, as is shown in Fig. 13 for the extended reionization 
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Figure 13: Rcionization shifts Cj^ to smaller angular scales by AZ ~ 20 around the region 
I ~ 50. For illustration Kr — 0.084 and r = 0.37 have been taken. 

Figure 14: Cf^ by the baryon isocurvature mode is positive around / ~ 100, whereas that 
by RGWs is negative there. Only when a very small r = 0.001 is taken, is the amplitude 
of Cf^ by the isocurvature modes comparable to that by RGWs. The Cf^ by isocurvature 
mode is the numerical result generated using CAMB [16] with q;_i = 0.015. . 

models with Kr — 0.084. Around the relevant region of I ~ 50, the reionization shifts the curve 
of Cf^ to smaller angular scales by an amount of AZ ~ 20, in comparison with the non-reionized 
Cf^ . This amount is much larger than that caused by NFS. Moreover, the shift AZ increases 
with the optical depth Kr- This significant effect has to be incorporated into the zero multipole 
analysis before one can make an extraction of RGWs from the total Cf^. 

In this procedure, besides disentangling the adiabatic (constant entropy) modes that are 
dominant in the scaler perturbations, one need consider the isocurvature modes possibly existing 
in the cosmological plasma [75], which can contribute to C^^^ [76]. In particular, the isocurvature 
modes contribute positively (correlation) to the cross spectrum Cf^ in the range around / ~ 100, 
in contrast to the adiabatic modes, which contribute negatively (anti-correlation) there. The 
observed Cf^ from WMAP has shown the anti-correlation, and a very stringent constraint has 
been found on the isocurvature contribution with the isocurvature/adiabatic ratio a_i < 0.015 
at 95%CL [5]. It is interesting to compare the contributions from RGWs and isocurvature 
perturbations to Cf^. The comparison is very sensitive to the ratio a_i and the tensor/scalar 
ratio r. Taking the upper limit q;_i — 0.015 constrained from WMAP-5, and using the CAMB 
Online Tool [16] results for isocurvature modes of the plasma components of baryon, CDM, and 
neutrino, one finds that when r = 0.37 is taken, the amplitude of Cf^ generated by RGWs is 
about two orders greater than that of the isocurvature modes. So in this case the isocurvature 
can be neglected. Only when a much smaller ratio r = 0.001 is taken, is the contribution by 
the isocurvature modes comparable to that by RGWs. This is demonstrated with r = 0.001 
and r = 0.01 in Fig. 14, in which the numerical Cf^ contributed by the baryon isocurvature 
perturbation has been produced from CAMB [16] with q;_i = 0.015. 

6. So far our analytic formulation for reionization can only distinguish two different extended 
models by comparing their Kr, rjr, and Arjr- The damping factor Dr{k) in Eq.(46) as a fitting 
formula could be used to specify other fine details of two reionization models. Obviously, with a 
fixed b, a larger c leads to lower bumps of Cf^ and C^^, as shown in Fig. 15 for the 2;-linear 



24 



Figure 15: The damping factor Dr{k) in Eq.(46) depends on the parameters c and b. The 
plot is made for the 2;-hnear model. 



model. On the other hand, with a fixed c, a larger b will yield a slightly higher reionization bumps, 
while leaving the primary peaks almost intact. This property can be inferred as the following. 
For the reionization bump around / ~ 5, the contribution is mainly from k ~ l/ivo — Vr) 
according to Eq.(73) , so Dr{k) oc e''^^^'^''^/^''"^'^'''''^^^ . In the reionization models considered in 
this paper, the combination Arirl/{rio — r]r) ~ 0.5 < 1, so a larger parameter b leads to a larger 
Dr and higher bumps. For the primary peaks with / > 100, Dr{k) is so small that the term 

2 

AlD?:(k) h(r]r) to CP^ and CP^ is practically negligible, leaving the primary peaks 

fe~i/(r)0-T?r) 

intact under a variation of b in Dr{k). 

6. Summary 

We have presented the approximate, analytical formulation of the reionized CMB spectra 
C/^^ generated by RGWs. Even though its approximate nature implies its application as a 
complement to the numerical codes, it does improve our understanding CMB and efficiently 
promote the analysis of various effects that reionization brings upon C^^. 

The reionization around 2; ~ 11 is studied by three simple homogeneous models, i.e., a sudden 
reionization, two extended reionizations with ionization fraction Xf,{r]) oc 77 and Xg{rj) oc z. The 
key parameter is Kr, the optical depth from the present back up to the start of reionization. 
Given a value of Kr in each model, the visibility function K-(^) follows, which is approximately 
fitted by Gaussian type functions. This procedure is similar to the treatment of decoupling in 
our previous study. 

Then the time integrations for polarization mode Pi and temperature anisotropies mode ai 
are carried out approximately, and the resulting analytic expressions consist of contributions by 
RWGs h{ria) and at the decoupling, and by h{rir) and h{rir) at the reionization as well. 

It is found that, while h{r]d) and h{rid) generate the primary peaks at / > 100, h{T]r) produces 
bumps for Cf and C[^^ at / 5, and h{r]r) enhances Cf^ and there. The analytic ^ 
qualitatively agree with those by the numerical computing, such as CAMB. 

As a merit of our analytic approach, the dependence of C/^^ upon the optical depth Hr are 
explicitly given, in terms of the coefficients ai{K,r) and a2{K,r) for the polarization Pi{r]o), and 
of the coefficients e~'^'' and (1 — e"'*'^) for the temperature anisotropies ai{r]o). It is found that 
ai{K,r) and a2{Kr) vary with Hr more quickly than e"'^'' and (1 — e"'^'^), respectively. Therefore, 
the polarization Pi is more sensitive to than the temperature anisotropies ai does. A larger 
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K,. gives higher a2{i^r) and lower ai(K,.), i.e., yielding higher bumps and lower primary peaks in 
Cf^ and Cj^^ . Thus there is a degeneracy of k,. with the normalization of the initial amplitude 
A of RGWs. Besides, Kr also has a weak degeneracy with the spectral index (3inf of RGWs since 
a larger jSinj enhances the primary peaks on small scales. The analytical Cf^ and Cj^^ also 
suggest possible ways to break these two kinds of degeneracies. 

Besides k,., our formulation also demonstrates the effects of the reionization time 77^ and the 
reionization duration Ar^^- For a fixed Kr, the height of bump is proportional to Ajy^, and the 
location I of bump depends on 77^ such a way I ~ k{riQ — rjr) that a later reionization (larger 
rir) yields a bump at larger angular scales (smaller I). 

Given a fixed set of parameters Hr, rjr, and Arjr, the 77-linear and z-linear models yield similar 
bumps in Cf'^ and C^^ . Thus our analytical formulation is unable to rediscover the reionization 
history Xe{r]) from C/^^. 

These analytical results tell that studies of reionization by means of CMB temperature 
anisotropics and polarization not only requires sufficient observational data, but also need detailed 
studies of the reionization process itself and more realistic modeling. 

The reionization process also significantly affects the possible detections of RGWs via the 
observations of C^-^ . In particular, it is found that the reionization causes a shift of the zero 
multipole Iq of the cross spectrum by a substantial amount A/ ~ 20, which is also Kr- 
dependent. The effect of reionization need be properly included, before one can apply the zero 
multipole method to extract the traces of RGWs from the observed Cf^. 

ACKNOWLEDGMENT: T.Y Xia's work has been partially supported by Graduate Student 
Research Funding from USTC. Y.Zhang's research work has been supported by the CNSF 
No. 10773009, SRFDP, and CAS. We thank Dr. Zhao for interesting discussions. 

References 

[1] D. N. Spergel, et al, Astrophys.J.Suppl., 148, 175 (2003). 

[2] H. V. Peiris, et al, Astrophys. J. Suppl., 148, 213 (2003). 

[3] D. N. Spergel, et al, Astrophys. J. Suppl., 170, 377 (2007). 

[4] L. Page, et al, Astrophys. J. Suppl., 170, 335 (2007). 

[5] E. Komatsu, et al, arXiv:0803.0547. 

[6] M.R. Nolta, et al, arXiv: 0803.0593. 



26 



[7] G. Hinshaw, et al, arViv:0803.0732. 

[8] J. Dunkley, et al, arXiv:0803.0586. 

[9] M. Zaidarriaga, Phys.Rev. D55, 1822 (1997). 

[10] M. Zaidarriaga, D. Spergel, U. Seljak, Astrophys. J. 488, 1 (1997). 

[11] M. Zaidarriaga and U. Seljak, Phys.Rev. D 55, 1830 (1997). 

[12] A. Venkatesan, Astrophys. J. 537, 55 (2000). 

[13] A. Venkatesan, Astrophys.J. 572, 15 (2002). 

[14] K.L. Ng and K.W. Ng, Astrophys.J. 456, 413 (1996). 

[15] U. Seljak and M. Zaidarriaga, Astrophys.J. 469, 437 (1996). 

[16] A. Lewis, A. Challinor and A. Lasenby, Astrophys.J. 538, 473 (2000). 
The CAMB Online Tool is located at 
http:/ /lambda. gsfc.nasa.gov/toolbox/tb_camb_form.cfm 

[17] S. Sasaki, Prog. Theor. Phys. 76, 1036 (1986). 

[18] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rep. 215, 203 (1992). 

[19] W. Hu and N. Sugiyama, Astrophys. J. 444, 489 (1995). 

[20] M.M. Basko and A.G. Polnarev, Mon. Not. R. Astron. Soc. 191, 207 (1980). 

[21] A. Polnarev, Sov. Astron. 29, 607 (1985). 

[22] R.Crittenden, D. Coulson, and N.G. Turok, Phys.Rev. D 52, 5402 (1995). 

[23] R. K. Sachs and Q. M. Wolfe, Astrophys.J. 147, 73 (1967). 

[24] B. Keating, P. Timbie, A. Polnarev, and J. Steinberger, Astrophys.J. 495, 580 (1998). 

[25] B. Keating, A. Polnarev, N. Miller, and D. Baskaran, Int.J. Mod. Phys. A 21, 2459 (2006). 

[26] U. Seljak and M. Zaidarriaga, Phys.Rev. Lett. 78, 2054 (1997). 

[27] M. Kamionkowski, A. Kosowsky and A. Stebbins, Phys. Rev. D 55, 7368 (1997). 

[28] A. Kosowsky, Annal Phys. 246, 49, (1996). 

[29] M. Zaidarriaga and D. Harari, Phys. Rev. D 52, 3276 (1995). 

[30] L. Grishchuk, Phys. Rev. D 48, 3513 (1993); Phys. Rev. Lett. 70, 2371 (1993); arXiv: 
astro-ph 0707.3319. 

[31] W. Zhao, D. Baskaran, and L.P. Grishchuk, Phys. Rev D 79, 023002 (2009). 

[32] J. Pritchard and M. Kamionkowski, Annals. Phys. 318, 2 (2005). 

[33] W.Zhao and Y. Zhang, Phys.Rev. D 74, 083006 (2006). 



27 



[34] T.Y. Xia and Y. Zhang, Phys.Rev. D 78, 123005, (2008). 

[35] S. Weinberg, Phys.Rev. D69, 023503 (2004). 

[36] D.A. Dicus and W.W.Repko, Phys.Rev. D 72, 088302 (2005). 

[37] Y. Watanabe and E. Komatsu, Phys.Rev.D 73, 123515 (2006). 

[38] H. X. Miao and Y. Zhang, Phys. Rev. D 75, 104009 (2007). 

[39] Y. Zhang, et. al, Class. Quant. Grav. 22, 1383 (2005); Y. Zhang, et. al, Class. Quant. 
Grav., 23, 3783 (2006); Y. Zhang , et. al, Int.JMPD Vol 17, 1105 (2008). 

[40] S. Wang, Y. Zhang, T. Y. Xia, and H. X. Miao, Phys.Rev. D 77, 104016 (2008). 

[41] D. J. Schwarz, Mod.Phys.Lett.A13, 2771 (1998). 

[42] L. Grishchuk, Lect. Notes Phys., 526, 167 (2001). 

[43] A.R. Liddle and D.H. Lyth, Phys.Lett. B291, 391 (1992). 

[44] L. Verde, et al., Astrophys.J.Suppl. 148, 195 (2003). 

[45] D. Baskaran, L.P. Grishchuk, and A.G. Polnarev, Phys. Rev. D 74, 083008 (2006). 

[46] Y. Cen, Astrophys. J. 591, 12 (2003). 

[47] T. Giannantonio and R. Crittenden, Mon.Not.R.Astron.Soc. 381, 819 (2007). 

[48] L.P.L. Colombo, et al., A^A 435, 413 (2005). 

[49] P.J.E. Peebles, Astrophys.J. 153, 1 (1968). 

[50] B.J.T. Jones and R.F.G. Wyse, A&iA, 149, 144, (1985). 

[51] A.J. Benson, A. Nusser, N. Sugiyanna, C.G. Lacey, Mon.Not.R.Astron.Soc. 320, 153 (2001). 

[52] G.C. Liu, et al, Astrophys.J. 561, 504 (2001). 

[53] M. Santos, et al., Astrophys.J. 598, 756 (2003). 

[54] J. Zhang, L. Hui, and Z. Hainnan, Mon.Not.R.Astron.Soc. 375, 324 (2007). 

[55] M. Mortonson and W. Hu, Astrophys.J. 657, 1 (2007). 

[56] P. Naselsky, L.Y. Chiang, Mon.Not.R.Astron.Soc. 347, 795 (2004). 

[57] S. Furlanetto and A. Loeb, Astrophys.J. 634, 1 (2005). 

[58] L. Hui and Z. Haiman, Astrophys.J. 596, 9 (2003). 

[59] X.H. Fan, C.L. Carilli, and B. Keating, Ann. Rev. Astron. Astrophys. 44, 415 (2006); X. 
Fan, et. al, Astron. J. 132, 117 (2006). 

[60] P.J.E. Peebles, Principles of physical cosmology, ( Princeton University Press, Princeton, 
1993). 

[61] R. Trotta and S. Hansen, Phys.Rev.D 69, 023509 (2004). 

28 



[62] L.M. Griffiths, D. Barbosa, and A.R. Liddle, Mon.Not.R.Astron.Soc. 308, 854 (1999). 

[63] L.P.L. Colombo, JCAP 0403, 003 (2004). 

[64] P. Mukherjee and A.R. Liddle, Mon.Not.R.Astron.Soc. 389, 231 (2008) 

[65] L.M. Griffiths and A. R. Liddle, Mon.Not.R.Astron.Soc. 324, 769 (2001). 

[66] M. Kaplinghat, et al., Astrophys. J. 583, 24 (2003). 

[67] L.P.L. Colombo and E. Pierpaoli, New Astron. 14, 269 (2009). 

[68] G. Jungman , M. Kaminokowski, A, Kosowsky, D.N. Spergel, Phys.Rev. D 54, 1332 (1996). 

[69] E.J. Eisenstein, W. Hu, M. Tegmark, Astrophys.J. 518, 2 (1999). 

[70] L.P.L. Colombo and S.A. Bonometto, New Astron. 8, 313 (2003). 

[71] M.J. Mortonson and W. Hu, Phys.Rev. D 77, 043506 (2008). 

[72] D. Baumann, et al, arXiv: astro-ph 0811.3919; D. Baumann, et al, arXiv: astro-ph 
0811.3911; M. Zaidarriaga, et al, arXiv: astro-ph 0811.3918. 

[73] A.G. Polnarev, N.J. Miller, and B.C. Keating, Mon. Not. R. Astron. Soc, 386, 1053 (2008). 

[74] N.J. Miller, B.C. Keating, and A.G. Polnarev, arXiv: astro-phy 0710.3651. 

[75] M. Bucher, K. Moodley, N. Turok, Phys.Rev. D 62, 083508 (2000); 
M. Bucher, K Moodley, N Turok, Phys.Rev. Lett. 87, 191301 (2001) 

[76] K. Enqvist, H. Kurki-Suonio, J. Valiviita, Phys.Rev. D 62:103003, (2000); 
K. Enqvist, H. Kurki-Suonio, Phys.Rev. D 61, 043002 (2000). 



29 




0.05 0.06 0.07 0.08 0.09 0.10 0.11 




11 



11 



1.0- 



CO. 

T3 
c 
CO 

d~ 

_c 

in 

H — ' 

c 

0) 

"o 

0) 

o 
O 



0.6- 



0.4- 



0.2- 



0.0- 



\ 






















\ 

\ 

v 








e 


-K 

r 














s 


a (i 


r' 




















\ 

N... 




















> 
























✓ ■ 


>» 
















/ 


a2(i 










> 










/ 

/ 


















f 

/ 

/ 








1 ■ 


-e" 


> 










1 






















'■ 1 

0.0 


1 

0.1 


1 

0.2 


1 

0.3 


1 

0.4 


1 

0.5 



K 



o 



z-linear 
T|-linear 



1 

i A 



^ -i.a 



o 



o 



o 



